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Abstract 

In this paper, we consider a simplified Erickscn-Lcslic model for the nematic liquid crystal 
flow. The evolution system consists of the Navicr-Stokes equations coupled with a convective 
Ginzburg-Landau type equation for the averaged molecular orientation. We suppose that 
the Navier-Stokes equations are characterized by a no-slip boundary condition and a time- 
dependent external force g(t), while the equation for the molecular director is subject to 
a time-dependent Dirichlet boundary condition h(i). We show that, in 2D, each global 
weak solution converges to a single stationary state when h(£) and g(t) converge to a time- 
independent boundary datum hoc and 0, respectively. Estimates on the convergence rate 
are also obtained. In the 3D case, we prove that global weak solutions are eventually strong 
so that results similar to the 2D case can be proven. We also show the existence of global 
strong solutions, provided that either the viscosity is large enough or the initial datum is 
close to a given equilibrium. 

Keywords: Nematic liquid crystal flow, non- autonomous Navier-Stokes equations, time- 
dependent Dirichlet boundary condition, long-time behavior, Lojasiewicz-Simon inequality. 
AMS Subject Classification: 35B40, 35Q35, 76A15, 76D05. 

1 Introduction 

We consider the following hydro dynamical model for the flow of nematic liquid crystals 

v t + v- Vv-i/Av + Vvr = -AV • (Vd © Vd) + g(t), (1.1) 
V-v = 0, (1.2) 
d f + vVd = r/(Ad-f(d)), (1.3) 

in x M + , where Q, C M n (n = 2, 3) is a bounded domain with sufficiently smooth boundary T, 
v = (vi, ...,v n ) tr is the velocity field of the flow and d = (di, d n ) tr represents the averaged 
macroscopic/continuum molecular orientations in R™ (n = 2,3). ir is a scalar function represent- 
ing the pressure (including both the hydrostatic and the induced elastic part from the orientation 
field). The external volume force is represented by g. The positive constants v, A and rj stand for 
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viscosity, the competition between kinetic energy and potential energy, and macroscopic elastic 
relaxation time (Deborah number) for the molecular orientation field. Vd Vd denotes the 
n x n matrix whose (i,j)-th. entry is given by Vjd • Vjd, for 1 < i,j < n. We assume that 
f(d) = V ( j-f 1 (d) for some smooth bounded function F : R n — > R. In particular, one uses the 
Ginzburg-Landau approximation f(d) = ^(|d| 2 — l)d to relax the nonlinear constraint |d| = 1 
on molecule length (cf. [TTlfTB] ). 

System (|l.ip - (|1.3p was firstly proposed in [IB] as a simplified approximate system of the 
original Ericksen-Leslie model for the nematic liquid crystal flows (cf. [5|ll5j ) . Well-posedness of 
the autonomous version of system (|l.ip - (|1.3p (namely, with g = 0, no-slip boundary condition 
for v, and time-independent Dirichlet boundary condition for d) has been analyzed in [18] (see 
also [6j[T2l[T9] and, for different boundary conditions, [22]). For numerical approximation we 
refer to [21 [ I23 [ |2"3|. Problem (|l.ip - (|1.3p has also been investigated on a Riemannian manifold 
in [26], where the existence of a global attractor in the 2D case was proven. As far as the long- 
time behavior of the single trajectory is concerned, in [18], a natural question on the uniqueness 
of asymptotic limit for global solutions (to the autonomous system) was raised. This question 
was answered in [3T], where it is proven that each trajectory converges to a single steady state 
(cf. [25] for some generalization). The proof is based on a suitable Lojasiewicz-Simon type 
inequality (see [27], cf. also [TO] and references cited therein). 

The technically more challenging case of time-dependent Dirichlet boundary conditions for 
d has been recently analyzed in [TJ[3llU[8]. For instance, under proper assumptions on the time- 
dependent boundary condition and assuming that g = 0, the existence of global weak solution, 
the existence of global regular solution for viscosity coefficient big enough, and the weak/strong 
uniqueness were obtained in [JJ. Regularity criteria for solutions in the 3D case can be found 
in [8]. Besides, the presence of a time-dependent external force is allowed in [1] and existence of 
global and exponential attractors is proven in the 2D case. In this paper, we want to extend the 
results of [3T] to the non-autonomous case treated in [T]. Thus we consider system (|l.ip - (|1.3|) 
subject to the boundary conditions 

v(x, t) = 0, d(x, t) = h(x, t), (x, t) E T x R + , (1.4) 

and the initial conditions 

v| t= o=v (a;) with V • v = 0, d| t=0 = d (x), xGfi. (1.5) 

In the 2D case, we prove that each weak/strong solution converges to a single stationary state 
when h(t) and g(t) converge to a time-independent boundary datum hoc and 0, respectively. 
In the 3D case, we first show the eventual regularity of global weak solutions, and the existence 
of global strong solutions provided that either the viscosity is large enough or the initial datum 
is close to a given equilibrium. Then an analogous result on the long-time behavior as in 2D is 
also obtained. In both cases, we provide an estimate on the convergence rate. 

Before ending this section, we state some key ingredients of the present paper. System 
(|1.1|) - (|1.5|) is non-autonomous due to the time-dependent boundary data h and external force 
g. This brings some additional difficulties into our subsequent proofs. First, in order to obtain 
the energy inequalities that play crucial roles in the proof of well-posedness as well as in the 
long-time behavior of global solutions (cf. Lemmas 12.21 12.51 12.61 15. ip . we have to introduce 
proper lifting functions (cf. (|2,7p and (|2.2ip below). The idea was first used in [3, 4J, but the 
lifting function introduced in this paper is different from the one in [3]. This is due to the 
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fact that we need some specific energy inequalities which not only yield uniform estimates of 
the solutions, but also provide estimates of the convergence rate (cf. Section 4). The second 
issue regards the application of the Lojasiewicz-Simon approach (cf. |27j ) which has been shown 
to be very useful in the study of long-time behavior of global solutions to nonlinear evolution 
equations (cf., for instance, [U|[TU | [TH |[50H52] and references therein). In particular, convergent 
results related to various evolution equations with asymptotically autonomous source terms were 
established, e.g., in [21171111)1111] . However, our current case is more complicated than the previous 
cases, because the Lojasiewicz-Simon inequality involves d which is subject to a time-dependent 
boundary datum. To overcome this difficulty, we derive an extended Lojasiewicz-Simon type 
inequality for vector functions with arbitrary nonhomogeneous Dirichlet boundary data, which 
is associated with the lifted energy (cf. Corollary 13. ip . This generalizes the results in [i~0ll31~] and 
may have its own interest. Third, in the 3D case, we also apply the Lojasiewicz-Simon approach 
to prove the existence of global strong solutions provided that the initial datum is close to a local 
minimizer of the elastic energy and the non-autonomous terms are properly small perturbations 
of their asymptotic limits (cf. Section 5). Then we further discuss the stability of these energy 
minimizers. This extends the previous results in |181l31j for the autonomous system, where the 
initial datum was required to be sufficiently close to a global energy minimizer. 

The remaining part of the paper is organized as follows. The next section is devoted to 
report some existence and uniqueness results and basic a priori estimates for the solution. The 
extended Lojasiewicz-Simon inequality we need is discussed in Section [3j In Section [J] we show 
the convergence of each global weak/strong solution to a single steady state and provide uniform 
estimates on the convergence rate in 2D. Results in 3D are presented in Section [SJ In particular, 
we study the eventual regularity of global weak solutions as well as the well-posedness when the 
initial data are close to local minimizers of the elastic energy. Long-time convergence of global 
solutions and stability of such minimizers are also proved. In the final Section [6j some useful 
properties of the lifting functions are reported. 

2 Preliminaries: well-posedness and a priori estimates 

Without loss of generality, from now on we set A = rj = 1. Let us introduce the function spaces 
we shall work with. As usual, L P (Q) and W k ' p (Q) stand for the Lebesgue and the Sobolev spaces 
of real valued functions, with the convention that H k (Q) = W The spaces of vector- valued 

functions are denoted by bold letters, correspondingly. Without any further specification, || • || 
stands for the norm in L 2 (ft) or L 2 (f2). This norm is induced by the scalar inner product 
(u, v ) = uvdx, where for vector valued functions the product uv is replaced by the Euclidean 
inner product u • v. We set, as usual, 

H = V L2(n) , V = V H(1)(Q) , where V = {v G C °°(O, R n ) : V • v = 0} . 

For any Banach space B, we denote its dual space by B*. In particular, we denote the dual 
space of Hj(fi) by H -1 (fi). 

In the following text, we will use the regularity result for Stokes problem (see, e.g., [29]) 

Lemma 2.1. For the Stokes operator S : D(S) = V n H 2 (S1) ->■ H defined by 

Su= -Au + Vvr G H, Vu G D(S), 
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it holds 

II u IIh 2 + IKIIifiyR < C||^ u ll) ^ u G 
for some positive constant C only depending on and n. 

We begin to report the existence of a weak solution (see [TJ Corollary 1.1]). 

Proposition 2.1. For any given T > 0, assume 

gGL 2 (0,T;V*), (2.1) 
h G L 2 (0,T;H§(r)), (2.2) 

h t GL 2 (0,T;H-^(r)) (2.3) 
|h| R n < 1, a.e. on T x [0,T], (2.4) 
d |r = h| i=0 . (2.5) 

Then for any (vo, do) G H x H 1 (J7) with |do |k« < 1 almost everywhere in Q, problem (|l.ip - (|1.5p 
admits a weak solution (v, d) suc/i 

v G L°°(0, T; H) n L 2 (0, T; V), 

d G L oo (0,T;H 1 (O))nL 2 (0,r;H 2 (^)), 

|d| K n < 1, a.e. on SI x [0,T]. (2.6) 
In the case n = 2, we also have the following results on uniqueness. 

Proposition 2.2. (cf. [TJ Theorem 1.4]) Le£ n = 2 and /ei i/ie assumptions of Proposition [Ql 
hold. Then the weak solution (v, d) to problem (|l.ip - (|1.5p giuen oy Proposition \2.1\ is unique. 
Moreover, we have (v,d) G C([0,T];H x H^O)). 

In order to obtain proper energy inequalities, we recall that suitable lifting functions were 
introduced in [3pE] to overcome the technical difficulties related to the time-dependent boundary 
datum for d. The first lifting function d E = &e{x, t) is of elliptic type (cf. [1]): 



\-Ad E = 0, inSlxR+ 
\d s = h, onrxR+. 

In particular, we define the lifting function d^o for the initial datum: 

(-Ad E0 = O, inn, 
)d so = d , on T. 



(2.7) 



(2. 



Set now 

d = d-d£. (2.9) 
Then system (|l.ip - (|1.5p can be rewritten into the following form: 

v t +v- Vv-i/Av + Vvr = -Ad-Vd + g(t), (2.10) 

V-v = 0, (2.11) 

dt + v-Vd = Ad - f (d) - d t d E (t), (2.12) 
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with homogeneous Dirichlet boundary conditions and initial conditions 

v = 0, d = 0, onrxK + , (2.13) 
v| t=0 = v , d|t =0 = d - d E0 , in ft. (2.14) 

Note that we have used the identity V • (Vd Vd) = ±V (|Vd| 2 ) + Ad • Vd to absorb the 
gradient term into pressure (cf. [IS])- 
Let us introduce the lifted energy 

£(t) = ^||v(t)|| 2 + \\\Vd(t)f + ^F(d(t))dx, t > 0. (2.15) 

Then we can derive the basic energy inequality for system (|l.l|) - (|1.5p . 

Lemma 2.2. Let the assumptions of Proposition \2.1\ be satisfied for all T > 0. Then, any weak 
solution which is smooth enough satisfies the following inequality for t > 

j2(t) + ^||Vv|| 2 + |][Ad - f(d)|| 2 < \\\8 t d E \\ 2 + C\\d t d E \\ + C\\g\&* := r(t), (2.16) 
where C is a positive constant independent of v and d. 

Proof. Multiplying (|2.1l)p and (|2.12[) by v and — Ad + f(d), respectively, integrating over ft and 
adding the results together, we get 

d (1 ll9 1 
di \r 

= (^d s ,Ad) + (g,v). (2.17) 



|v|| 2 + i||Vd|| 2 + j F{d)d^\ +H|Vv|| 2 + ||Ad- f(d)|| 2 



In above, we have used the facts (v • Vv, v) = (VP, v) = (v • Vd, f(d)) = due to the 
impressibility condition V • v = 0. By the Poincare inequality ||v|| < Cp||Vv|| and (|2.6p . the 
right-hand side of (|2.17p can be estimated as follows 

|(d t d E ,Ad) + (g,v)| 

< | (d t d E , Ad - f (d)) | + | (d t d E , f (d)) | + | (g, v) | 

< ||Ad - f(d)||||a t d E || + ||f(d)||||^d E || + ||v|| v ||g||v 

< ^l|Vv|| 2 + |]|Ad - f(d)|| 2 + ^||^d B || 2 + C||^d E || + C||g||^. 

The proof is complete. □ 

Remark 2.1. We fix the calculations in Q Lemma 2] in which the term (dtd E , Ad) is missing. 
Though it does not affect the proof of existence, it does have influence on the long-time behavior 
of global solutions (especially the convergence rate). 

Let us now introduce the following (Banach) spaces of translation bounded functions 
L q tb (0, +oo; X) :={h G L q loc ([0, +oo);X) : 

rt+l 

ll h ll^(o,+oo ; x) : = sup jf ||h(r)||^dr < +oo| 



where X is a (real) Banach space and q £ [1, +oo) is given. 

From the basic energy inequality (|2,16p . through a suitable Galerkin approximation scheme, 
one can derive uniform-in-time estimates for any weak solution (the proof is a minor modification 
of [U Lemma 1.2, Remark 1.1]). 
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Lemma 2.3. Let the assumptions of Proposition \2.1\ hold for all T > 0. In addition, suppose 
that 



gGL 2 (0,+oo;V*), (2.18) 
hGL? 6 (0,+oo;Hl(r)), (2.19) 
h t G L 2 (0,+oo;H~5(r)) n +oo;H - 3(T)). (2.20) 

Then a weak solution (v,d) to problem (|l.ip - (|1.5p given by Proposition \2. 1\ is a global solution 
on [0, +oo) and fulfills the following uniform bounds 

||v(i)||<C, ||d(t)|| H i < C, Vt>0, 

/Vl|Vv(r)|| 2 + ||(Ad - f(d))(r)|| 2 )dr < C, V t > 0. 
J o 

Here C is a positive constant depending on ||vo||, ((doHn 1 ; l|g||L 2 (o,+oo ; v*); ll n ll 2 3 , 

||h+|| „ i and \\hA\ i 

1 n L 2 (0,+oo;H"2(r)) 11 "Li^+oojH"^ (T)) 

Next, we introduce the lifting function dp = dp(x,t) of parabolic type, which satisfies 

'<9 t dp-Ad P = 0, in fi x M+, 
dp = h, on T x R+, (2.21) 
dp(0) = dp , in 0. 

The motivation of introducing the parabolic lifting function dp is that we now have, by definition, 
A(d — dp) — f(d)|r = 0. This fact is crucial when we use integration by parts to derive some 
higher-order differential inequalities of system (|1.1|) - (|1.5|) (cf. [HE]). We note that dp in (|2.21|) 
is different from the one introduced in [5] as they have different initial values. Both choices are 
valid for the proof of existence result, but the current definition of dp is necessary for the study 
of long-time behavior. Denote 

d = d - dp. 

System (|l.ip - (|1.5p can now be rewritten into the following form: 

v 4 + V'Vv-i/Av + VP = -Ad-Vd + g(t), (2.22) 
V-v = 0, (2.23) 
d t + vVd = Ad-f(d), (2.24) 

with homogeneous Dirichlet boundary conditions and initial conditions 

v = 0, d = 0, on r x R+, (2.25) 
v|f= = v , d|t= = d - dp , in CI. (2.26) 

In the sequel, we shall frequently use the following lemma (cf. [1]) 
Lemma 2.4. The following equivalence between norms hold 

||v||ni ~ ||Vv||, ||d|| H i f» ||Vd||, in Hg(fi), 

||v|| H 2 ps ||Av||, |d||jj2 ps |Ad||, mHj(S])nH 2 (n), 
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||d|| H 3 rj ||V(Ad)|| + ||Ad||, in Hj(n) n H 3 (ft). 

7/d and dp are functions which are smooth enough and |d|jjn < 1, |dp|]R« < 1, then we have 

||Ad|| < ||Ad P || + ||Ad-f(d)|| + C, 
||VAd|| < ||VAdp|| + ||V(Ad- f(d))|| + C||Vd||. 

where C is a positive constant independent of d and dp. 

Let us introduce the quantity 

^ P (t) = ||Vv(t)|| 2 + ||Ad(t) - f(d(t))|| 2 , t > 0. 

Lemma 2.5. Let n = 2 and Zei i/ie assumptions of Lemma \2.3\ hold. If the weak solution (v, d) 
is smooth enough then it satisfies the following inequality 

^p(t) < C(A 2 P (t) + A P (t) + (2.27) 

where 

fli(t) = ||^dp(t)|| 4 + ||^dp(t)|| 2 + ||VAdp(t)|| 2 + ||g(t)|| 2 . (2.28) 

HereC is a positive constant depending on v , \\vq\\, ||do Hh 1 ? llg|lL 2 (o,+oo ; v*); ll^ 1 *ll i 2( 0+ H _ ^(r))' 
||h+|| , i and llhll „ 3 

11 ^'-^(O.+ooiH-^r)) 11 M L2 b (0,+oo;H2(r)) 

Proof. Taking the time derivative of Ap(t), we obtain by a direct calculation that 

±±A P (t) + (v\\Svf + (Ad- f(d)\\ 2 ) 

= -(Sv, v Vv) + (Sv, g) - (Sv, Ad • Vd) - (V(v • Vd), V(Ad - f (d))) 

-(f(d)d t ,Ad-f(d)) 

5 

:= (2.29) 
i=i 

To get this identity we have used the fact that Ad — f (d)|r = as well as (Sv, v t ) = (— Av, v 4 ). 
It is not difficult to see that 

l-Til < ||5v||||v|| L 4||Vv|| L 4 < C||5v||(||Vv||5|| v ||5)(||Av||^||Vv||5) 

< C||Sv||l||Vv|| < e||Sv|| 2 + C||Vv|| 2 , 
\h\ < e||Sv|| 2 + C||g|| 2 . 

For ^3, we have 

\I 3 \ = |(Sv,(Ad-f(d))-Vd) + (Sv,f(d)-Vd) + (Sv,d t dp-Vd)| 

< ||5v||||Vd|| L 4||Ad-f(d)|| L 4 + ||Sv||||Vd|| L oo||^dp|| 

< e||Sv|| 2 + C||Vd|| 2 4 ||Ad - f(d)|| 2 4 + C||Vd|||oo ||cW|| 2 - 

On account of Lemma 12.31 we infer from the Sobolev embedding theorems that 

||Vd|| 2 4 < C||Ad||||Vd|| + C||Vd|| 2 < C||Ad|| + C||$d P || + C 
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< C||Ad-f(d)|| +C\\d t d P \\ +c, 
||Vd||£cc < C||VAd||||Vd|| + C||Vd|| 2 

< C||VAd P ||+C(l + ||V(Ad-f(d))||), 
||Ad-f(d)|| 2 4 < C||V(Ad-f(d))||||Ad-f(d)||. 

Using the above estimates, we obtain the estimates for I3 and I4: 

\h\ < e||5v|| 2 + C||V(Ad - f(d))||||Ad - f(d)||(||Ad - f(d)|| + ||d t d P || + 1) 
+C||d t dp|| 2 (||VAd P || + ||V(Ad - f(d))|| + 1) 
< e||5v|| 2 + e||V(Ad - f(d))|| 2 + C||Ad - f(d)|| 4 + C|| Ad - f(d)|| 2 
+C||^dp|| 2 (||a t d P || 2 + ||VAdp|| + 1), 

\h\ < ||V(Ad-f(d))||(||Vv|| L 4||Vd|| L 4 + ||v|| L ^||d|| H 2) 

< e||V(Ad - f(d))|| 2 + C||Vv|| 2 4 ||Vd|| 2 4 + CUvH^dlAdH 2 + 1) 

< e||V(Ad - f(d))|| 2 + C||Av||||Vv||(||Ad - f(d)|| + ||9 f d P || + 1) 
+C||Av||||v||(||Ad - f(d)|| 2 + HdidpH 2 + 1) 

< e||Sv|| 2 + e||V(Ad - f(d))|| 2 + C||Vv|| 4 + C\\Ad - f(d)|| 4 
+C||Vv|| 2 + C||Ad - f(d)|| 2 + ||d t d P || 4 . 

We now observe that 

Is = -(f'(d)d t , Ad - f (d)) - (f'(d)^dp, Ad - f (d)) := I 5a + h h . (2.30) 
Recalling (|2.24p . we have 

|/ 5o | = |(f'(d)(v ■ V)d, Ad - f (d)) - (f'(d)(Ad - f(d)), Ad - f (d))| 

< ||f / (d)|| L ^(||v|| L 4||Vd|| L 4||Ad - f(d)|| + ||Ad - f(d)|| 2 ) 

< C||Vv||||v||||Vd|| 2 4 +C7||Ad- f(d)|| 2 

< C||Vv|| 2 + C||Ad - f(d)|| 2 + C\\d t d P \\\ 

\hb\ < ||f / (d)|| LO o||9 t dp||||Ad-f(d)|| <C||Ad-f(d)|| 2 + C||9 i dp|| 2 . (2.31) 

Finally, collecting the above estimates and taking e sufficiently small, we deduce that 

j t A P (t) + (u\\Sv\\ 2 + ||V(Ad - f(d))|| 2 ) 
< C(A 2 P (t) + A P (t)) + C\\d t dp\\ 2 (\\d t d P \\ 2 + ||VAdp|| + 1) + ||g|| 2 , 

which easily implies the inequality (|2.27p . □ 

Taking advantage of Lemmas 12.51 16.11 and 16.21 one can deduce the following results on the 
regularity of weak solutions as well as the existence of strong solutions to system (|l.ip - (|1.5p in 
2D. 

Theorem 2.1. Let n = 2 and let the assumptions of Proposition \2.1\ hold for all T > 0. In 

addition, suppose that 

g G L 2 (0,+oo;H), (2.32) 
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h G L t 2 fe (0, +oo;H 2(r)), (2.33) 
hi e L 2 (0,+oo;H5(r)) n L^O, +00; H _ 2(r)). (2.34) 

(i) System fjl.Xj) — f|1.5|) admits a unique global weak solution (v, d) satisfying 

||v(t)|| v < (^(l + t" 1 ), ||d(t)|| Ha < (^(l + t- 1 ), Vt>0, 
j\\W(.r)\\ 2 H2 + ||d(r)||^ 8 )dr < C(l + Or, t e [S,T], 

where C is a positive constant depending on v , ||vo||, ||d.o Hh 1 j Hg|lL 2 (o,+oo;H)> l|h|| 2 s , 

('mJ 7f(vo,do) £ V x H 2 (0), i/ien problem (|1.1|) - ([1.5|) admits a unique global strong solution 
(v, d) satisfying 

\\v(t)\\ v <C, ||d(i)|| H 2 < C, Vi>0, (2.35) 
/(l|v(r)||^ + ||d(r)||^)dr<C7r, te[0,T], (2.36) 

JO 

whereC isapositive constant depending onv, ||vo||v> ||do||H 2 > l|g||L 2 (o,+oo;H) > ||h 
llh II 1 llh II 1 

11 *"L 2 (0,+oo;Hi (r))' 11 i "i 1 (0,+oo;H _ 2(r))' 



L 2 6 (0,+oc;H?(r))' 



Remark 2.2. Lemma \2.3\ and Theorem \2.1\ still hold when g and hj are translation bounded 
with respect to time (see \T$). 

Next, we consider the 3D case. Instead of Lemma 12.5} we have the following higher-order 
energy inequality 

Lemma 2.6. Let n = 3 and let the assumptions of Lemma \2.'J\ hold. If a weak solution (v, d) 
is smooth enough then it satisfies the following inequality 

±A P (t) + (u- Cl A P (t)) l|Sv|| 2 + - ^A P (t)) l|V(Ad - f(d))|| 2 
< C(l + u~ 2 )(A P (t) + R 2 (t)), t>0, (2.37) 

where Ap(t) = Ap(t) + 1 and 

R 2 (t) = \\d t d P (t)\\ 2 + \\d t d P (t)\\ 6 + ||Va t d P (t)|| 2 + ||g(t)|| 2 . (2.38) 

Here c±,C2,C are positive constants that may depend on ||vo||, Udolln 1 an d on l|g||L 2 (o,+oo:V*); 
llhll „ , 3 , , llhtll „ 1, , ||ht|| , , 1, „, but they are independent of v. 

11 "z4(o,+oo;H2(r))' 11 tM L 2 (0,+oo;H-^(r))' 11 1 11 L 1 (0,+oo;H _ 2 (r)) ' " 1 J 

Proof. We estimate the right-hand side of (|2,29p by using the 3D version of Sobolev embedding 
theorems. We have 

\h\ < ||-SV||||v|| L 6||Vv|| L 3 < C||5v||||Vv||(||Av||5||Vv||5) 
< C||SV||i||Vv||5 < ^||Vv||i||Sv|| 2 + C||Vv|| 2 , 

h < a 5v 2 + - gf. 

8 v 
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Recalling that ||d|| H i < C (cf. Lemma l2.3|) . from the Sobolev embedding theorems as well as 
Agmon's inequality in dimension three, we infer 

||Vd|| L3 < C||Ad||l||Vd||5+C'||Vd||<C'(||Ad-f(d)|| + ||^d F ||)i+C', 
||Vd|| L6 < C||Ad|| + C||Vd|| <C\\Ad-f{d)\\+C\\d t d P \\+C, 

||Vd|| L « < C||Vd||| 1 ||Vd||| ? <C'(||VAd||5||Ad||3 + ||Ad||+l) 

< C(||V(Ad-f(d))||3||||Ad-f(d)||3 + ||Ad - f(d)||*||VAd P ||3 
+||V(Ad - f(d))||3||Ad P ||3 + ||VAdp||5||Adp||3 + ||V(Ad - f(d))||s 
+ ||VAd P ||5 + || Ad - f(d)|| + ||Adp|| + 1), 

||Ad-f(d)|| L 3 < C||V(Ad-f(d))||3||Ad-f(d)||i 

Thus we have 

1^1 < ||SV||||Vd||i,a||Ad-f(d)|| L 3 + ||5v||||Vd|| L cc||^dp|| 

< C||5v||(||Ad - f(d)|| + ||^dp|| + l)||V(Ad - f(d))||5||Ad - f(d)||3 
+C\\Sv\\\\d t d P \\ (||V(Ad-f(d))||3||Ad-f(d)||3 + ||Ad-f(d)||3||V5tdp||5 

+ ||V(Ad - f(d))||5 \\dtdp\\$ + ||V^dp||3||^d P ||3 + ||V(Ad - f(d))||3 

+||Va t d P ||3 + ||Ad - f(d)|| + ||5 t d P || + l) 

< ^ + i||Ad - f(d)|| 2 ) ||5v|| 2 + i||v(Ad - f(d))|| 2 

+C(1 + v- 2 )(\\Ad - f(d)|| 2 + ||9 t dpf + ||d t d P || 2 + ||VAd P || 2 ), 

|J 4 | < ||V(Ad-f(d))||(||Vv|| L 3||Vd|| L 6 + ||v|| L oo||d|| H 2) 

< C||V(Ad-f(d))||||Vv||3||Av||3(||Ad-f(d)|| + ||^dp|| + 1) 

- (i + |n vv ii a ) n Sv ii 2 + Q + ^rii A5 - f ^n 2 ) n y ( Aa - f ™i 2 

+C||Ad - f(d)|| 2 + C{1 + Ol|Vv|| 2 + C||5 4 dp|| 4 , 

\ha\ < ||f / (d)|| LO o(||v|| L 6||Vd|| L 3||Ad-f(d)|| + ||Ad-f(d)|| 2 ) 

< C||Vv|| 2 ||Vd|| 2 : +C||Ad-f(d)|| 2 

< C||Av||||v||(||Ad - f(d)|| + \\d t d P \\ + 1) + C||Ad - f(d)|| 2 

< ^ll^v|| 2 + C(l + 0(||Vv|| 2 + || Ad - f(d)|| 2 + ll^dpH 2 ). 

o 

We observe that I$b can be estimated as in (|2.3ip . Then, collecting all the estimates of Ij, we 
have 

^p(0+(^-||Vv||i-||Vv|| 2 ) l|5v|| 2 
+ ^l__L||Ad-f(d)|| 2 ) ||V(Ad-f(d))|| 2 
< C(l + z,- 2 )(^lp(t) + H^dpH 2 + H^dpf + ||VAdp|| 2 + ||g|| 2 ). 
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As a result, there exist constants c\,C2 > independent of v such that the following inequality 
holds 



^ P (t) + (1/ - Cl ^ P (t))||5v|| 2 + M - ^p(t)) HV(Ad - f(d„ 

< C(l + i/- 2 )(^4p(i) + ||a t dp|| 2 + H^dpf + ||VAdp|| 2 + ||g|| 2 ), 

which implies (|2.37|) . □ 

On account of Lemma [2.6| one can deduce that system (jl.ip - (jl.5p admits at least one global 
strong solution, provided that the viscosity is large enough (see jH Theorem 7] for the case g = 0, 
cf. also [T8j[31] for the autonomous case). We just report a result under weaker assumptions 
than in [1] and omit the detailed proof. 

Theorem 2.2. Let n = 3 and assume that (|2.32p - (|2.34p and (|2.4p are satisfied. For any 
(vo,do) G V x H 2 (fi) satisfying (|2.5p and |do|]R3 < 1, there exists a Vq > 0, depending on 

||(Vo,do)|| V xH» «^l|g||^(0,+oc;H)> H h H i?i)( o, +00; H§ (F)) ' l|h * 'W+oo;^ (F)) ' h^+oo^^ (F)) ' 

suc/i that, for any v > uq, problem (|l.ip - (|1.5p admits a global strong solution (v,d) which sat- 
isfies the same uniform estimates as in the 2D case (cf. (|2.35|) and (|2.36p ). 

Remark 2.3. When n = 3, the weak-strong uniqueness result obtained in J^], Theorem 7] still 
holds in our case. Thus, the global strong solution (v,d) obtained in Theorem \2.2\ is unique. 

3 Extended Lojasiewicz— Simon type inequality 

For all d G M := {<j> G H 1 (il) : 0|r = hoc}, where G H2(T) is given, we consider the 
functional 

E(d) = l\\Vd\\ 2 + / F(d)dx. (3.1) 



2 

It is straightforward to verify that 

Lemma 3.1. If i]> G H 1 ^) is a u>eaA; solution to the elliptic problem 

[ 0|r = hoc, 

i/ieri ^ is a critical point of the functional E(d) in J\f. Conversely, if tp is a critical point of the 
functional E(d) in J\f , then if) is a weak solution to problem (|3.2p . 

3 

Remark 3.1. If is more regular, then ip is more regular. For instance, if G H2(r), 
then ip G H 2 (S1). 

Lemma 3.2. Suppose that ip is a critical point of E(d) inN '. Then there exist constants (3\ > 0, 
6 G (0, |) depending on ip such that, for any w G N that satisfies ||w — V'IIh 1 < Pl> there holds 

|| - Aw + f(w)|| H -i > \E(w) - E(ip)\ l - Q . (3.3) 



Remark 3.2. The above lemma can be viewed as an extended version of Simon's result \21^ for 
scalar function under I? -norm. We can refer to 11 (X Chapter 2, Theorem 5.2], in which the 
vector case subject to homogeneous Dirichlet boundary condition was considered. We observe that 
the result can be easily proved by modifying the argument in U0\j using a simple transformation 
(cf. also 122 Remark 2.1]). 
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The Lojasiewicz-Simon type inequality ()3.3p only applies to proper perturbations of the 
critical point of energy E in the set J\f and it is not enough for our evolutionary problem (jl.ip - 
(jl.5p . whose boundary datum is time-dependent (not necessary in M). In order to overcome this 
difficulty, we prove the following extended result that also involves the perturbation of boundary: 

Theorem 3.1. Suppose that ip is a critical point of E(d) in N . Then there exists a constant 
ft £ (0, 1) depending on ip such that, for any d £ H 1 (r2) satisfying ||d — ip\\ni < /3, there holds 



C \\d\r — hoc 1 1 i , 

ii °°iiH2(r) 



+ ||d| r - h^ll 1 ^ ) + || - Ad + f (d)||H-i > \E(d) - Ei^l 1 ' 



(3.4) 



where 8 £ (0, |) is the same constant as in Lemma WJ^ while C is a positive constant depending 
on ip. 



Proof. For any d £ H 1 (r2), we have that Ad £ H Then we consider the elliptic boundary 

value problem 

f Aw = Ad, 
| w| r = hoo. 

It easily follows from the elliptic regularity theory (cf. e.g., [551 Proposition 5.1.7]) that 



(3.5) 



|w - d|| H i < C||d| r - hoo|l H i (r) , 



(3.6) 



which implies 



w 



| H i < ||w - d|| H i + ||d 
< C||d-V|| H i- 



| H i <C||d| r - hoc|l H i (r) + ||d 



H 1 



(3.7) 



Let f3i be the constant in Lemma 13.21 We infer from the above inequality that if /3 £ (0, 1) is 
chosen sufficiently small, then we have ||w — V'IIh 1 < Pi- As a consequence of Lemma 13,21 we 
have 

|| - Aw + f(w)|| H -i > \E(w) - Eitf;)^ 6 . (3.8) 



Aw + f (w)||h-i < II - Ad + f (d)||H-i + C||f(d) - f(w)|| L e (n) 



On the other hand, by the definition of w, we can see that 

^(w)-^)! 1 ^ < II- 

< || - Ad + f(d)|| H -i +C||d- w|| H i 

< || - Ad + f(d)|| H -i + C||d| r - hoo|l H i (r) - 

We deduce from £ (0, \) that 

\E(d) - E^)] 1 ' 6 < \E(w) - E^l 1 - 9 + \E(d) - ^(w)l 1 - 

and 



(3.9) 



(3.10) 



\E(d) - E(w) 



< 



1-6 



iVdl 



iVwl 



1 1—0 




+ 


/ 

Jn 



(F(d) - F(w))dx 



< C(||d|| H i, ||w|| H i)||d - wf-, 6 < C||d| r - h oo||^ (r) , 



where in (|3,lip we use the facts that ||d|| H i < HV'Hh 1 + P an d ||w|| H i < HV'Hh 1 + Pi- 
Hence, combining ()3.9p — (j3. lip , we deduce (|3.4p . 



(3.11) 



□ 
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Since the basic energy inequality (|2.16p (cf. Lemma l2.2p is only valid for the lifted energy £ 
(|2.15p . in order to apply the Lojasiewicz-Simon approach to our problem, we need to consider 
the following auxiliary functional corresponding to energy E (cf. (|3.ip ): 



E(d) = -||Vd|| 2 + I F(d)dx, Vdetfffi), 
2 Jn 



(3.12) 



where 



d = d — d#, 

and d^ is the elliptic lifting function satisfying the following elliptic problem (cf. ([2 

j-Ad E = 0, 
\d£; = d|r, xer. 

Then we have 



(3.13) 



Corollary 3.1. Suppose that ip is a critical point of E(d) in N . Then there exist constants 
ft € (0, 1) and 9 G (0, ^) depending on ip such that, for any d G H 1 (f2) satisfying \\d — i/'||h 1 < P> 
there holds 



C\\d\ T - hooll 1 ^ + || - Ad + f(d)|| H -i > |#(d) - E^)\ l ~ 
HJ(r) 

where C is a positive constant depending on ip and . 
Proof. From the definition of E(d), we set, for ip G A/", 

%,) = I||V^|| 2 + f F(iP)dx, 



(3.14) 



(3.15) 



where ip = ip — ipE and ipE satisfies 



-AipE = 0, x G 0, 

V>£ = hoc, xer. 



(3.16) 



A direct calculation yields that 

E(d) = 

E{ip) = 



£(d) + -||vd £ || 2 



E(i/j) + -\\Vip E \\ 2 



Vd : Vd E dx, 
Vip : S/tpEdx, 



where we used the notation A : B = Y27j=i ^ijBij . Theorem 13.11 implies that there exist 
constants ft G (0, 1) and 9 G (0, |), such that for any d G H 1 (il) satisfying ||d — ip\\ni < ft, 
holds. Next, we proceed to estimate the quantity |-E(d) — E^ip)] 1-0 



\E(d)-EW\^ d 
< lEW-E^O+f^ 



V(d E -ip E ) :V(d E + ip E )dx 



+ 



(Vd : Vd E - Vip : Vip E )dx 



Ji + J 2 + h 



(3.17) 



13 



The estimate for J\ follows from (|3.4p . Since ||d — V'IIh 1 < /? < 1, then ||d|| H i < H^Hh 1 +1- For 
J2, we infer from the elliptic estimate (cf. [28} Proposition 5.1.7]) that 

h < C||V(d £ -^)|| 1 - e ||V(d i3 + ^)|| 1 ^ 



< C||d| r -h c 



H2(r) 



r 1 + hoo 1 
11 n H ^ (r) 11 °°iiH2(r) 



< C||d| r -hooll 1 ^ (lldl 

H2(r) 

< C||d| r -hooll 1- / ■ 

H2(r) 



Hi(n) + l|hoo|| H i (r) 



Recalling the function w introduced in ()3.5p . we estimate J3 as follows 



< 



[V(d - w) : Vd E + Vw : V(d E - ip E ) + V(w - ip) : Vip E ] dx 



Vw : V(d E - i) E \ 



+ 





l-e 




w) : Vd E dx 


+ 


J 






Jn 




l-e 





V(w — ip) : S/if) E dx 



'■— ha + hb + Re- 
using (|3.6p and (|3.7p and the fact ||d 



H 



1 < j3, we observe that 



ha < ||^(ci - w)!! 1 -^!!^^!! 1 - < c||ci| r - h — lli I "i (r) N d lrll^ (;r) 

< C||d|r-h TC ||XJd|Ko) < C H d lr-hoo|lX,' 



HJ(r) 



H2(r) 



(3.18) 



(3.19) 



(3.20) 



hb < ||Vw|| 1 ^||V(d i; -^)|| 1 - e <C7(||V||Hi+C/3) 1 - 9 ||d|r-h 0O || 1 - 1 9 

H5 (r) 

< C||d| r - hooll 1 - . 

H2(r) 



(3.21) 



For he using integration by parts and noticing that Aip E = 0, (w — ip)\r = 0, we obtain 

/ V(w - ij>) : Vtp E dx = - / (w - i/>) ■ A-ip E dx + / (w - V0|r ■ d n ip E dS = 0, (3.22) 
Jn Jn Jr 

where n is the unit outer normal to the boundary T. Thus (|3.22[) implies that 

he = 0. (3.23) 

Finally, since 1 — 9 G (0, 1), we have lldlp— hoJI 1 < Clldlr — hooH 1- ] 61 . Summing up, we can 

H2(r) HJ(r) 

conclude from ([33D, (l3TT)l - (l3T23]l . and Ad = Ad that (137L4]) holds. The proof is complete. □ 

Remark 3.3. If 9 G (0, ~) is such that (pUl"]) holds, then, for all 9' G (0,9) and any d G H^fJ) 

satisfying ||d — V'IIh 1 < P> we s ^// have 



C (\\d\ v - hooll 1 -, 6 ' + || - Ad + f(d)|| H -0 > \E[d) - E(n^ d ', 



(3.24) 
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where C is a (properly adjusted) positive constant depending on ip and h^. To see this, we first 

i _oi i-e 7 1-9' i-e' 

notice that, since 2 > i— ^ > 1, for any a,b>0, it holds (a + b) ^ < 2 (a ^ + b i-» ). TTien it 

follows from (|3.14|) t/iat 

i^(d) - %)r- e ' = (i£(d) - s wi 1 - 9 ) u 

1-9' 

< f ||d| r - hooll 1 " +||-Ad + f(d)|| H -i 

< cflldlr-hooll 1 - 6 ' +||-Ad + f(d)|| H -i 

V" H^(r) 




4 Long-time behavior in 2D 

In this section, we focus on the case n = 2. In order to study the long-time behavior of global 
solutions to problem (|l.ljl - (|1.5jl . we need some decay conditions on the time-dependent external 
force g and boundary data h, namely, 

(HI) f t + °° ||h t (r)|| H x (r) dr < (7(1 + t)" 1 ^; 
(H2) f+°° ||h t (r)|| 2 H , (r) dr < (7(1 + t)~^- 

(H3) !F\WT)t*r<C(l + t)-^ 
(H4) ||g(t)|| 2 <C(l + t)^; 
(H5) ||hi(t)|| L2(r) <C(l + t)^; 

for all t > 0. Here C and 7 are given positive constants. We also note that (H4) entails (H3). 

Since in the 2D case weak solutions become strong for positive times (cf. Theorem 12. ip . 
we can confine ourselves to consider strong solutions. We recall that, for any given global 
strong solution (v, d), we have the uniform estimate (|2.35p . It follows that the w-limit set of 
the corresponding initial datum (vo,do) is non-empty. Namely, for any unbounded increasing 
sequence {t n }^Li, there are functions Voo G V and doo G H 2 (f2) such that, up to a subsequence 
{tj}?° =1 C {t n }£Li, we have 

lim ||v(tj) - Vqo 1 1 = 0, lim ||d(tj) - doo|| H i = 0. (4.1) 

j— >+oo J— H-00 

Next, we characterize the structure of the w-limit set. In order to do that, we first recall a 
technical lemma (see [321 Lemma 6.2.1]) 

Lemma 4.1. Let T be given with < T < +00. Suppose that y and h are nonnegative con- 
tinuous functions defined on [0, T] and satisfy the following conditions: ■£ < c\y 2 + C2 + h, 
with y(t)dt < C3, Jq h(t)dt < C4, where Ci (i = 1,2,3,4) are given nonnegative constants. 
Then for any p G (0, T), the following estimates holds: y(t + p) < + C2P + c<^J e clC3 , for all 
t G [0,T — p]. Furthermore, ifT = +00, then lim y(t) = 0. 

t— >+oo 
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Proposition 4.1. Let the assumptions of Theorem \2. 1\ hold. Then the uj -limit set uj(vq, do) is 
a subset of 

S = {(0, u) : u G n H 2 (ft) siic/i i/iai - Au + f (u) = 0, in ft}. 

Moreover, we have 

lim ||v(t)|| v = 0, (4.2) 

t— >+oo 

lim ||-Ad(t) + f(d(t))||=0. (4.3) 

t— >+oo 

Proof It follows from Lemma E2 that / + °° ||Vv(t)|| 2 + ||Ad(i) - f(d(t))\\ 2 dt < +oo, which 
together with the definition of Ap and (|6.3p yields 

r+oo r+oo 

/ A P (t)dt< / (||Vv(t)|| 2 + 2||(Ad-f(d))(t)|| 2 + 2||a f d P (t)|| 2 )(it<+oo. (4.4) 
Jo Jo 

Using Lemma 12.51 and Lemma 14.11 we can see that 

lim A P (t) = 0, 

which implies lim t ^ +00 ||Vv(t)|| = 0. Hence, for any (vo^d^) G o;(vo,do), we have Vqo = 0. 
On the other hand, by definition of Ap, (|4.4|) also yields that 

lim ||-Ad(t)+f(d(*))||=0. (4.5) 

t— >+oo 

From Lemma 16.21 we have lim ||c?tdp(t)|| = (cf. (16.50 ). As a result, it follows from the 

t— s-+oo 

inequality 

0< ||- Ad(t) + f(d(t))|| < ||-Ad(i) + f(d(t))|| + ||d t d P (t)||, i>0 (4.6) 

that (|4.3p holds. Concerning the limit function doo, we infer from (|2.35p that doo G H 2 (ft) and 
(|4.ip holds. We now check the boundary condition for doo. Since G i 1 (0, +oo; H~2 (r)), h(t) 
strongly converges to a certain function hoc G H~2 (r) as time goes to infinity with a controlled 
rate, namely, 



Hh^-hooll^i^ < / ||h t (r)|| H _ i(r) dr-).0, ast^+oo. (4.7) 



On the other hand, we infer from (|233"j) and ([235]) that h G L°°(0, +oo; Ha (r)). Consequently, 

3 3 

hoc G H2 (r) and h weakly-* converges to h^o in L°°(0, +oo; H2 (r)). By interpolation, we have 
lim \\h(t) — h 00 || L 2c r - ) = 0. Thus, from the asymptotic behavior of the boundary datum h, we 

t— >+oo y ' 

have for any j G N, 

< ||doo|r — hoo|| L 2( r ) < ||doo|r — h(*?')llL 2 (r) + ||h(*i) ~~ n oo||L 2 (r) 

< CUdoo - d(^)|| H l + \\Htj) - Ml^T)- 

Hence, letting j — > +00 in the above inequality, we deduce from (|4.ip and (|4.7p that doo|r = hoo- 
For any z G Hj(ft) and j G N, we have 



/ (-Adoo + f (d^)) • zdx 
Ju 
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< f (-Adoo + Ad(i,)) ■ zdx + /" (f (<!«,) - f (d(tj))) • zda; 

+ / (-Ad(%) + f (dfo))) ■ zdx 

< ||V(d(tj) " doo)||||Vz|| + (C||d(tj) - doollni + II - Adfo) + f(d(^ 
Passing to the limit as j — > +00, we get 

f (-Ad oo + f(d oo ))-zdx = 0. 

Summing up, we can see that doo G W n H 2 (fi) solves (|3.2p . The proof is complete. 
We can also prove the convergence of the lifted energy. 



□ 



Proposition 4.2. Let the assumptions of Theorem \2. 1\ hold. Then the lifted energy functional 
£ defined by (|2.15p is constant on the uj-limit set cj(vo,do). Namely, there exists a constant <?eo 
such that E{d OQ ) = £00, for all (0, doo) with d ro G N C\ H 2 (f2). Moreover, we have 



lim £(t) = £ 

t->+oo 



(4. 



Proof. From the previous argument, we know that for arbitrary (0, doc?), (0, dj^/) G w(vo,do) 
there exist unbounded increasing sequences {t^}'jL 1 and {tj}'^ =1 such that (|4.ip holds. As a 
result, we have 

. lim £(tf) = £(d«), lim f (tf) = E(d£>). 
On the other hand, it follows from the basic energy inequality (|2,16p that for any t' > t" > 0, 



,(2> 



\£{t')-£(t")\< I r(t)dt->0, ast',t" 
Jt" 



+00. 



Then by 



|E(d«) - E(dW)\ < \£{tf) - £{tf)\ + \E{S£) - £{tf)\ + \£{tf) - B(dg>)|, 

letting j — > +00, we can see that E(d^) = E(d^). Namely, £ is a constant (denoted by £oo) 
on the w-limit set cj(vo, do). Moreover, for any t > there exist tj < tj + i such that t G [tj, tj+i] 
and |£(t) - ^l < |£(t) - £ (tj)\ + - fool which yields (|Q]> . □ 



4.1 Convergence to equilibrium 

Theorem 4.1. Lei i/ie assumptions of Theorem \2.1\ hold. If, in addition, we assume (H1)-(H3), 
then any strong solution (v(t),d(£)) convergence to an equilibrium (0, doo) strongly in VxH 2 (r2) 
as t goes to +00. 

Proof. On account of (|4.2[) we only need to prove that d{t) converges to doo as t — > +00 given 
by (|4.ip . Below we adapt the idea in [2,7J to achieve our goal. Indeed, observe that we can find 
an integer jo such that for all j > jo, \\d(tj) — doo|| H i < |, where j3 G (0, 1) is the constant 
given in Corollary 13.11 (depending on doo). Consequently, we define 

s(tj) = sup{r > tj : ||d(r) - doollH 1 < P}- 
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Since d £ C([0, +00); H 1 (0)), we can see that s(tj) > tj for any j > jo. By Lemma 12.21 and 
Proposition 14.21 we have 

\S(t) - S(doo)| > - min{^, 1} y P 2 (r)dr - jf r(r)dr, 

where 

P(t) = ||Vv(t)|| + ||Ad(t)-f(d(t))||, 
and r is defined in (|2.16p such that, thanks to (H1)-(H3), we have 

/+00 
r(r)dT < C(l + fT 1_7 , Vt>0. 

Let the constant 9 be as in Corollary 13.11 (depending on doo)' Using Remark 13.3] we can choose 
9 f E (0, 9] such that 6 f also satisfies 

< " < WTTY m 

If 9 itself satisfies f|4.9[) . we just take 0' = 0. For any fixed tj with j > jo, we introduce the sets 

Kj = fo, %)), Kf> = [t € K : V(t) > (1 + £ )-(i-^')(i+t) | ) ^P) = ^ \ R (i) 
Consider the following functional on Kj 

§(t)=£(t)-E(d 00 ) + 2j r(r)dT, VteKj. 

It easily follows that 



lim $(t--) = 0. (4.10) 

j— S-+CO 



Next, we have 



4(|*(*)fsgn*(*)) = iU( t ) < Jl min{l/) l}\$(t)f-i-D 2 (t) < 0, (4.11) 

at at 4 

which implies that the functional |$(i)| e 'sgn<I>(i) is decreasing on Kj. Keeping in mind that 
9' < 9 and 2(1 — 9') > 1, we can apply Corollary 13.11 (cf. also Remark I3.3P to obtain that 

f r+oo \ 1-0' 

/-.N 2(1-0') _ / , +00 \l-e' 

< (ij l|v|| 2(1 " e,) + C||h(t) - hooll^jf + Cll - Ad + f (d)||H-i +C[J r(r)drj 

(r+oo \ 1-^' / />+oo \ 1-0' 

_/ l|ht(r)|| H i (r) drJ +Cfy r(r)drj 

< C||Vv|| + C\\ - Ad + f(d)|| +C(l + t)" (1 - 9 ' )(1+7) . (4.12) 



Thus, on -ft'j 1 ^, we have 
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which together with (|4.1ip yields that on 

-±(Mt)\ 6 'sgn$(t))>CV(t). (4.13) 

As a consequence, we have 

/ v(t)dt <-c I j t {m)f *mm)dt < cmt^f + ms^f) < +oc, (4.u) 

J J 

(2) 

where <&(s(tj)) = if s(tj) = +00. On the other hand, on K- , we have 

p poo /~i 

/ V{t)dt<C (l + t)- ( - 1 - e '^ 1+ ^dt = {l + t j f e ' +0 '-\ (4.15) 

Here, we notice that j9' + 0' - 7 < due to (|P]> . Then (gUj) and (|4"T5j) imply that 

/ V(t)dt = / i V{t)dt + / 2 V(t)dt < +00, 

■> 3 3 

for any j. On the other hand, it follows from (|2.35|) and (|2.12|) that 

||d t (*)|| < ||v • Vd|| + || Ad - f(d)|| < ||v|| L 4||Vd|| L 4 + ||Ad - f(d)|| < CV(t). (4.16) 
As a consequence, 

/ \\d t (t)\\dt < cm tj )f + Msit^f) + c{\ + t j y e ' +e '-\ (4.i7) 

Proposition 4.3. Let the assumptions of Theorem \2. 1\ hold. Then there exists an integer j\ > jo 
such that s(tj 1 ) = +00. Thus 

||d(t) - doollni < A yt>t h . 

Proof. The conclusion follows from a contradiction argument (cf. |13|). Suppose that for any 
j > jo we have s(tj) < +00. Then, by definition, we have 

||d(a(ti)) - doollni = P > 0. (4.18) 

Besides, it follows from (|4TT]l . (jlTOD and (OT) that 

||d(s(t i ))-d 00 || < ||d(a(t,-)) — d(tj)|| + ||d(t i )-d 0O || 

rs(tj) 

< / ||d t (*)||cft + \\d(tj) -doo|| -> 0, asj^+oo. 



Using uniform estimate (|2.35p and interpolation inequality, we obtain 

||d(s(ij)) - doollni < \\d(s(tj)) - d 00 || H2 ||d(s(tj)) - d^H -> 0, as j -> +00, 
which leads a contradiction with (|4,18p . The proof is complete. □ 
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Due to Proposition 14.31 we have = +00 for some j\ > jo. Arguing as above, we can 

prove 

r+00 

/ \\d t {t)\\dt < +00. 
Thus d(t) converges in L 2 and recalling (|4.ip . by compactness we conclude that 



lim ||d(t) - doollni = 0. (4.19) 

t— >+oo 



Finally, observe that 



|Ad(t) - Adooll = || - Ad(t) + f (d(*))|| + ||f (d(<)) - f(doo)[| 

< ||-Ad(t) + f(d(t))|| + C||d(t)-d 00 ||. (4.20) 



Then (03]) and (|4"T9|) entail that 



lim ||d(i) - doollna = 

t— >+oo 



and this finishes the proof. □ 
4.2 Convergence rate 

Theorem 4.2. Let the assumptions of Theorem \2. 1\ hold. If, in addition, we assume (H1)-(H2) 
and (H4)-(H5), then we have 

||v(f)|| + ||d(t) - d^llHi < C(l + ty^w, t > 0. 
Moreover, if (H2) and (H5) are replaced by, respectively, 
(H6) \\h t (t)\\ Hi{r) <C(l + t)-^; 

(El) HKtJ-hooll^^^CCl + t)- 1 -^ 
the following higher-order estimate holds 

||v(t)|| v + ||d(t) - dooHna < C(l + ty^r, t > 0. 

Proof. The proof consists of several steps. 

Step 1. L 2 -estimate of d — doo. This follows from an argument devised in [7]. For the 
readers' convenience, we sketch the proof here. From the previous argument, we only have to 
work on the time interval [tj 1 , +00). Denote 



Since 



/+00 
r(r)dT. 



j t Ht) < -- min{i/, l}V\t) - r(t) < 0, 



and lim <J>(t) = 0, we know that <J>(t) is decreasing and §>(t) > for t > t,- . 

First, if the boundary datum h and the external force g become time-independent in finite 
time, i.e., there exists time Tq such that for t > To, h = and g = 0. Then the problem 



20 



reduces to the autonomous system considered in [31] . Thus, below we just assume that either 
h or g does not become time-independent in finite time (namely, the system will always be 
non-autonomous). In this case, if there exists t* > tj 1 such that <J>(i*) = 0, then T>(t) = r(t) = 
for all t > t* and this is a contradiction since r(t) cannot identically vanish from any finite time 
on. Therefore, we can suppose 

$(t)>o, yt>t h . 

If the open set is bounded, then there exists t* > tj 1 such that [t*, +oo) C As a 

result, V(t) < (1 + t)^( 1 - e ')(i+7) an d by (ETTBj) . we have 

/+oo S~i 
W^)\\dr < _ 7g/ _ g/ + 7 (i + ty e ' +e '-\ v t > t*. 

Next, we treat the case when the open set is unbounded. There exists a countable 
family of disjoint open sets (a n ,b n ) such that = U^L 1 (a n , b n ). On , recalling (|4.12j) . 
we can see that on any (a n , b n ) C k£\ it holds 



As a result, for any t G (a n , 6 n ), 



< ^(«n) ie - 1 + C(l - 26K)(* - an)j l - 2e , (4.21) 

where by the definition of and (|4.12p we have 

$K) < CD(a n )^ + C(l + a n )~( 1+ ^ = C{1 + any 1 "''. 
Using the fact (1 + 7)(1 — 20') > 1 (cf. (|4.9p ). we can take n* £ N sufficiently large such that 

<S>(a n *f e '- 1 - C(l - 20'K* > a^ 1-2 ^ - C(l - 20'K* > 1. (4.22) 
Therefore, we infer 

$(*) < (^(l + t)"^ 7 , V t G (a n ,,oo)n^] i 1) . 
Similar to (|4. 13|) . we have (since <&{t) > 0) 

- j t m 9 ' > CV(t), V t G (a n *,oo) n 

Due to (|4.9p . it follows that — 7$' — 0' + 7 > tSoqt- Now for any t > a n *, we can conclude that 

r+00 r p 

\\d(t) - dooll < / ||d t (r)||dr = / ||d t (r)||dr + / ||d t (r)||dr 

J(i,oo)nitj i 1) J{t,oo)nKf> 

r r+00 

< C V{r)dT + C (l + r)-^- d '^ 1+ ^dT 



< c<s>(tf + c(i + ty e ' +e '-~< 

e' 

< C(l + i) i- 29 '. 
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Using (|2.35p . after properly adjusting the constant C, we have 

||d(t) -doo|| < C(l + t)~^ 7 , \ft>0. (4.23) 
Step 2. H x H 1 -estimate. It easily from the basic energy inequality (|2.16p that 

j t y{t) + ^l|Vv|| 2 + |||Ad- f(d)|| 2 < r (t), (4.24) 

where 

y(t) = \h{t)\\ 2 + \\\v(d(t) - doo)!! 2 + j[F(d){t) - F(d 00 ) - f(doc)(d(t) - d^dx. 

As in [31], using (|2.35p . we can show that 



f [F(d)(t) - ^(doo) - f(doo)(d(t) - d^dx 



< C\\d(t) -d 



|2 

oo I ■ 



Keeping in mind the definition of lifting functions, we have d — doo|r = so that 

HVCd-doo)!! < cHACd-doo)!! < || - Ad + f (d)[| + a||f (d) - f (doo)!! 

< C|| - Ad + f (d)|| + C||d(*) - dooll, 

HVCd-doo)!! < ||v(d-d 0O )|| + ||v(d f; -d 0O )|| 

/+oo 
||h t (r)|| H i (r) dr. 

Thus it follows that 

11 / f' +QC \ 2 

y(t) > -||v(t)|| 2 + -||V(d-d oo )|| 2 -C||d(0 -d^ll 2 -CU ||h t (r)|| H i (r) drJ(4.25) 

y(t) < C||Vv|| 2 + C|| - Ad + f(d)|| 2 + C||d(t)- dooll 2 . (4.26) 

Condition (JO) implies that < 7- Then we deduce from (|4T24l . (EL"23]) . (H4)-(H5) and 

Lemma 16.11 that 

j t y(t) + ay(t) < C(r(t) + \\d(t) - d^f) < C(l + 4)"^, (4.27) 
where a > is sufficiently small. The above inequality implies that 

26' 

y{t)<C{l + t)~^r, Vt>0. (4.28) 
Combining it with (|4.25p and recalling (HI), we get 

J ll h *WII H ^(r) dr ) 

28' 

< C(l + t) T ^ T , \ft>0. (4.29) 

Step 3. V x H 2 -estimate. Taking advantage of the stronger assumptions (H6)-(H7) and 
(|4.29p . we now get a higher-order estimate. Observe first that 

|| - Ad + f(d)|| < || - Ad + f(d)|| + ||Ad P || = || - Ad + f(d)|| + HfljdpH, 
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then we have 

y(t) < C||Vv|| 2 + C\\ - Ad + f(d)|| 2 + C||^d P || 2 + C\\d(t) - dooll 2 . 

It follows from ^TFf\ and fljjZZD that 
d 



dt 



z(t) + a 2 z(t) < C(l + ty^w + C(i?i(t) + \\dtdp 



where 

«(t) = j(t) + Mp(<), 

and a\ and «2 are sufficiently small positive constants. From the definition of R\, ()6.6 
fact ^7 < 2 + 2 7 (cf. gSD), we have 



26' 



i*i(t) + \\d t d P (t)\\ z < C(l + + c||VAd P (t)|r. 



Hence from (|4.30p we infer that 



z(t) < z(0)e- a2t + Ce~ a2t f e a21 

Jo 



< Ce- a2t + e~ a2t I e a2T 



+e~ a2t I e a2T 



C(l + t)~t=w + ||VAd P (r)|| 2 



dr 



C(l + r) + ||VAd P (r 

20' 



C(l + r )-T3w + ||VAdp(r)|| 2 

:= Ce- a2t + Z 1 {t) + Z 2 (t). 
It follows from ([63]) and (H6) that 



a 2 -1 / 2 



C(l + r)~^2F + ||VAdp(r) 



< Ce-^ t \t+ \l + T)- 2 - 2 " l dT)<C{l + t) 



Next, by (|6.7p and the fact 2 \ d i < 1 + 27, we deduce that 



+ ^ j ± e a2T dr + C £ \\VAd P (r)\\ 2 dT 



Z 2 (t) < Ce~ a2t [l + - 

28' . _ 26' 

< c(i+ty^r +c(i + t)~ 1 - 27 < cii + ty^ 7 . 

As a result, we obtain that 

26' 

z(t) < C(l + t) i-2«', V f > 0. 

In particular, we have 

26' 

A P (t) < C(l + t) T ^ w , t > 0, 
which together with ()4.20p and (|4.29p yields the following estimate 

||v(t)|ft + ||Ad(t) - Ad^H 2 < C(l + ty^, t > 0. 
Finally, using a standard elliptic estimate, we obtain (cf. (H7)) 

l|d(t) - d 00 || H 2 < C\\Ad(t) - Ad^ll + C\\h(t) - hoc || 3 „ < C(l + ty^w 

H2 (1 j 

for all t > and this finishes the proof. 



(4.30) 

(4.31) 
and the 

(4.32) 



(4.33) 



(4.34) 



(4.35) 

(4.36) 
(4.37) 



□ 
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5 Long-time behavior in 3D 

As in the classical Navier-Stokes case (see [14j). we can prove the eventual regularity of any 
global weak solution. Thus the convergence results can also be extended to the 3D case. Indeed, 
comparing with Lemma 12.61 we derive first an alternative higher-order energy inequality. 

Lemma 5.1. Let the assumptions of Proposition [Ol hold for all T > 0. Suppose, in addition, 
that (|2,32|) - (|2.34|) are satisfied. If a weak solution (v,d) is smooth enough then it fulfills the 
following inequality 

j t A P {t) + v\\Svf + l|V(Ad - f(d))|| 2 < CMW) + A P {t) + R 3 (t)), (5.1) 

where 

R 3 (t) = \\d t d P (t)\\ e + \\d t d P (t)\\ 2 + ||VAd P (t)|| 2 + ||g(t)|| 2 , (5.2) 
for allt > 0. Here C* is a positive constant that may depend on v, ||vo||, Udolln 1 ; IIsIIl 2 (o,+oo;V*)> 

" h "L t 2 6 (0,+oo;Hi(r))' " h *"L 2 (0,+oo;H-i(r)) ; " ^ "i 1 (0,+oo;H- h (r)) ' 

Proof. We reconsider the estimates in the proof of Lemma 12.61 Recalling (|2.29p and (|2,30D , 
thanks to the Young inequality, it is not difficult to obtain that 

|ii| < ||SV||||v|| L 6||Vv|| L 3 < C||5v||i||Vv||i <e||5v|| 2 + C||Vv|| 6 , 
\h\ < e||Sv|| 2 + C||g|| 2 , 

\h\ < e||^v|| 2 +e||V(Ad-f(d))|| 2 +C||Ad-f(d)|| 6 + C7||Ad-f(d)|| 2 

+C||d t d P || 6 + C||a t dp|| 2 + C||VAdp|| 2 , 
\h\ < e||Sv|| 2 + e||V(Ad - f(d))|| 2 + C||Ad - f(d)|| 6 + C||Vv|| 6 + C||Vv|| 2 + C||d f d P || 6 , 
\ha\ < e||5v|| 2 + C||Ad-f(d)|| 2 + C||Vv|| 2 + C||9 i dp|| 2 . 

In addition, 1^ can be exactly estimated as (|2.3ip . Collecting all the estimates, and taking e to 
be sufficiently small, we obtain our conclusion (|5.1|) . □ 

Then we prove the following sufficient condition. 

Proposition 5.1. Suppose that the assumptions of Proposition flTTl and ()2.32p - (|2.34p are satis- 
fied. In addition, assume that (vo,do) £ VxH 2 (il). If there exists a sufficiently small Eq G (0, 1] 
such that 

H|Vv(*)|| 2 + ||Ad(t) - i(d(t))\\ 2 )dt < e . (5.3) 

'o 

then problem (|l.ip - (|1.5p admits a unique global strong solution (v, d) in Q x (0, +oo), provided 
that llhdl i is small enouqh. 

Proof. For simplicity, we give a formal proof. To make it rigorous we should work within a 
proper approximation scheme (see, for instance, [US]). Let Li > (i = 1,2,3,4,5) be the 
constants such that 



/ 

Jo 



v || + ||d || H i < Li, (5.4) 

l h *lli, 2 (0 >+ oo;H^(r)) ^ L2 ' ^ 

< L 3 , (5.6) 



'Z&(0,+oo;H3(T)) 
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Ilhtll , i, < L 4 , (5.7) 

l|g||L 2 (0,+oo;V*) < L 5- (5-8) 

It follows from the basic energy inequality ()2. 16|) that 

£( t ) + \J^ (^llVvf + ||Ad- f(d)|| 2 ) <£(0) + J^ r(t)dt, Vt > 0. (5.9) 

Then, by definition of E and Lemma 16,14 we have 

l|v(t)|| + ||d(t)|| H i <Ci, Vt>0, (5.10) 
J ^||Vv|| 2 + ||Ad-f(d)|| 2 J^<C72, (5.11) 

where the constants C\, C2 depend on L±, L5 and 
Let K > be such that 

HlVvoH 2 + ||Ad(0) - /(d )|| 2 < K. (5.12) 
Keeping Lemma l5.1l in mind and arguing as in [18], we consider the following Cauchy problem 

^-Y(t) = C,(Y(t) 3 + Y(t)) + C*R 3 (t), Y(0) = max{l,u~ 1 }K>A P ( y 0). (5.13) 

We denote by / = [0,T max ) the (right) maximal interval for the existence of a (nonnegative) 
solution Y(t) so that lim Y(t) = +00. It easily follows from (|5.ip and the comparison principle 

that < Ap{t) < Y(t), for any t £ I. Consequently, Ap(t) is finite on I. We deduce from 
Lemma 16.21 that 

r+00 

Rs(t)dt < C 3 , 

where C3 is a constant depending on Q, ||g||L2( 0> +oo;H) an d ^2- Besides, we note that T max 
is determined by Y(0), C* and C3 such that T max = T max (Y(0), C*, C3) is increasing when 
y(0) > is decreasing. Taking to = \T max > 0, then it follows that Y(i) (as well as Ap(t)) is 
uniformly bounded on [0, to]. This easily implies the local existence of a unique strong solution 
to problem (|1.1|) - (|1.5|) (at least) on [0, to] (actually on [0,T max ), but we lose uniform estimates 
on such maximal interval). 

By Lemma 16.21 (cf . (|6.15p ) , we have 

sup \\A(d P (t) - d E (t))\\ 2 < c\\h t \\ 2 r (5.14) 

t>0 L 2 (0,+oo;H2(r)) 

where c is a constant that depends only on f2. Set now 

toK) T _______ |\ t2 K 



e = min<{ l,-^- } , L 6 = min<J 1,L 2 , — } . (5.15) 



From the assumption, there exists a small constant £0 < £ o such that ()5.3[) is satisfied. Therefore, 

to 
2 

H|Vv(^)|| 2 + ||Ad(t») - f(d(^))|| 2 < 2e i " 1 - 



we can find t* € [y,io] such that 
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Moreover, if we further assume 



then by ()5.14p we obtain 



|hj 2 i < L 6 , 

' "L 2 (0,+oo;H2(r)) - 



A P {U) < i/||Vv(t*)|| 2 + ||Ad(t*)-f(d(tA)|| 2 

< u\\Vv(Q\\ 2 + 2||Ad(t*) - f(d(t*))|| 2 + 2||A(d P (t*) - d E (t,))\\ 2 

< l /||Vv(^)|| 2 + 2||Ad(t,) -f(d(^))|| 2 + 2c||h t || 2 o 

L 2 (0,+oo;H2(r)) 

< 4e to 1 + 2cL 6 < K 

< max {1, i/ -1 } if = Y(0). 

Taking t* as the initial time for equation f|5. 13j) . we infer from the above argument that Ap(t) 
is uniformly bounded at least on [0, ^] C [0,t* + to]- Moreover, its bound only depends on Q, 
v, Li, ...,Lq, C* and to- Then by an iterative argument we can show that Ap(t) is uniformly 
bounded for all t > and this enable us to extend the local strong solution to the whole time 
interval [0, +oo). The proof is complete. □ 

A consequence of the above proposition is the eventual regularity of global weak solutions. 

Theorem 5.1. Suppose that the assumptions of Proposition [Ol and (|2.32p - (|2.34|) are satisfied. 
Let (v,d) be a global weak solution to (|1.1|) - (|1.5|) . Then there exists a large time T* G (0, +oo) 
such that (v,d) is a strong solution on (T*,+oo). 

Proof. Let L\, L2, L3, L4, L5 > be the constants as in the proof of Proposition 15.21 For a weak 
solution (v,d), we still have the uniform estimates (|5.10p and (|5.1ip . Considering the ODE 
problem (|5.13p . we can fix the constants e~o, Lq and to- Taking £0 = £q, we observe that there 
must exist a sufficiently large T\ > such that 

H|Vv|| 2 + || Ad- f(d)|| 2 ) dt < e , (5.16) 

Ti 

||Ad P (t)- Ad £ (t)|| < L 6 , Vt>[Ti,+oo), (5.17) 

where for the second inequality we have used Lemma f6.2f i) and the fact that dtdp(t) = Adp(t) — 
Ad E (t). Also, (pU6l) implies that there is T* G [T x ,Tx + 2t ] such that 

H|Vv(7;)|| 2 + ||Ad(T*) - f(d(T,))|| 2 < ^. (5.18) 

As a result, 

H|Vv(TA|| 2 + ||Ad(T*) - f(d(TA)|| 2 

< H|Vv(7;)|| 2 + 2||Ad(TA - f(d(TA)|| 2 + 2||A(d P (T*) - d s (T,))|| 2 

< ^°+2cL 6 

< K. 

Taking T* as the initial time, then we can apply Proposition 15 . II to conclude that problem (jl.ip 
(|1.5p admits a unique global strong solution (v,d). By the weak/strong uniqueness result [U 
Theorem 7], we see that (v, d) coincide with (v,d) on [T*,+oo). The proof is complete. □ 
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Thanks to the eventual regularity we can argue as in the previous section to prove 



Theorem 5.2. Suppose that the assumptions of Theorem I5.il hold. Then any global weak so- 
lution given by Proposition \2.1\ converges in V X H 2 (il) to a single equilibrium (0,doo) with 
estimates on the convergence rate similar to the 2D case, provided that g and h fulfill the cor- 
responding hypotheses (H1)-(H7) as in Theorems \4-l\ and\4-2\ 



Remark 5.1. We recall that there exists a (unique) global strong solution when the viscosity 
is large enough (cf. Theorem \2.2}) . Consequently, due to Lemma \2.6\. all the results proven in 
Section^ (i.e., Theorem \4-l\ and Theorem \4-2\ ) still hold with the same assumptions on the data. 
The related proofs just require some minor modifications. 

The existence of a global strong solution is also ensured (with no restrictions on viscosity) 
when the initial data are close to a given equilibrium and the time dependent boundary data 
satisfies suitable bounds. First, recall that the basic energy inequality (|2.16p implies (cf. (|5.9p ) 

ft r+oo 

/ M|Vv(t)|| 2 + ||Ad(t) - f(d(t))\\ 2 )dt < 2(f (0) - £{t)) + 2 / r(t)dt, 
Jo Jo 

and 

r(t)dt < C r ( ||h t || 2 x + llhJI i „ + ||g||?2m^-v*i I > ( 5 - 19 ) 

where C r is a universal constant. Then we can easily deduce from Proposition 15.11 that if the 
lifted energy stays sufficiently close to its initial state, then system f[l.ip — (|1.5|) admits a unique 
global strong solution (cf. [18] for the autonomous case) . 

Proposition 5.2. Assume (|2.32p - (|2.34p and (|2.4j) hold. Moreover, suppose that (vo,do) £ 
V x H 2 (f2) satisfying (|2.5p and |do|]R3 < 1. If there exists a sufficiently small Eq € (0,1] such 
that 

£(t) > £(Q) - e , Vt > 0, (5.20) 

where £ is the lifted energy defined by (|2.15p . then problem (jl.ip - (|1.5p admits a unique global 
stronq solution (v,d) mOx(0,+oo), provided that II h+ II „ i , ||hd| i and 

l|g||i 2 (o,+oo;V*) are small enough. 

We can prove the global existence of a strong solution that originates near a local minimizer 
of the lifted energy. For this purpose, we assume that for all t > (comparing with assumptions 
(HI), (H4), (H5)) 

(HP) /+~ ||h t (r)|| H x (r) dr < Ahil + t)-^; 

(HP) ||g(t)|| 2 <M 2 (l + t)- 2 -T; 
(H5') \\h t (t)\\ L 2 {r) <M 3 (l + t)- 1 ^. 

Here Mj, j = 1,2,3 and 7 are positive constants. 7 characterizes the decay rate of non- 
autonomous terms, while Mj control their magnitude. 
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Theorem 5.3. Assume (|2.32p - (|2.34|) and flZ3]) hold. Moreover, suppose that (v ,d ) <E V x 
H 2 (0) satisfying (|2.5p and |do|n3 < 1. Denote by d* E the unique solution to 

f-Ad*=0, xell, 

< S (5.21) 

[d^ = h^, x e r, 



and set 



(d) = i||V(d - d^)|| 2 + / F(d)dx, Vd g AA. 



Lei d* € A/"n H 2 (0) 6e a Zoca/ minimizer of £{d) in the sense that ${d) > ${d*) for all d £ N 
satisfying \\d — d*|| H i < 5. Suppose also that the initial data vo and do satisfy 

||v || v <l, ||d -d*|| H 2 < 1. (5.22) 

There exist positive constants o~i, o~2, M±, M2, M3, Lq, which may depend on the system coeffi- 
cients, on and on d* , such that if the initial data (vo,do) and h also fulfill 

llvnll < fi, lldn — d*|| w i < (Jo, ||ht|| 2 1 < Ln, 

II un _ i, n u 11 h _ z, 11 f ll L2(0)+oo . H i (r)) - u ' 

and (HI'), (H4'), (H5') hold with such Mj, j = 1,2,3, and 7 > 1, then problem (fTI] ) -(|T3 j) 
admits a unique global strong solution (v, d). 

Proof. Without loss of generality, we assume 5 £ (0,1]. In the subsequent proof, Q (z G N) 
stand for a positive constant which only depends on O, z/, 7 and d*. Since ()5.22|) holds, it is not 
difficult to see that the constants L\ and K in (|5.4|) and ()5.12|) depend on d* only. We just take 
L2 = L3 = L4 = L5 = 1 in (|5.6p for the sake of simplicity. Then we have the uniform estimate 
(cf. (foTTUD ) 

llvWII + l|d(*)ll H i < Ci, i>0. 

Arguing as in the proof of Proposition 15.11 we find that problem (jl.ip - (jl.5p admits a unique 
strong solution (at least) on [0, to], whose V x H 2 norm is uniformly bounded on [0,io] : 

||v(t)|| v + ||d(t)||H? < C 3 , i € [0,t ]. (5.23) 

Besides, we can also fix the constants £0 an d L§ (see (|5.15p ). Here, we just take £q = Eq and 
Lq = Lq. It follows from (|5.19p that 

f+00 



p +00 

/ r(t)dt < C r C s (Mi + M 2 + M|) < -p , 



provided that Mi,M2,M3 > are assumed to be properly small and satisfying 

Mi + M 2 + M 3 2 < £ ° 



4C r C ' 



where C s is a universal constant due to the Sobolev embedding. Hence, according to Propositions 
15.11 and 15.21 in order to prove the existence of global strong solution, we only have to verify that 

£(t)-S(0) > -y, Vt>0. (5.24) 
First, we notice that (recalling ([22]), (|23]) and (ETT2]) ) 
£(0)-£(t) < ^|vo|| 2 + £(d )-£(d(t)) 
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= ^llvoll 2 + i||V(do - d E0 )\\ 2 - i||V(d(t) - d E )f + J F(d ) - F(d(t))cfc 
< ^l|vo|| 2 + C4(||do-d(t)|| H i +||d B0 -d s || H i). (5.25) 
On the other hand, thanks to standard elliptic estimates, we have 

||dEo-dE|| H i < c||d |r -h(t)|| H i 

< c ||do|r - hoo|l H i (r) + c||hao - h(t)|| H i (r) 



< c||do-d*|| H i +cj ||h t (r)|| H i (r) dr 

< ca 2 + cMx, Vt>0. (5.26) 

Let 

a^minjl,^}, a 2 < ^ min{l, c~% M 1 <min{l,^}. 
Due to (|5.25p and (|5.26|) . in order to prove (|5.24|) . we only have to verify 

||d -d(i)|| H i<^, Vt>0. (5.27) 

Since d* G jVn H 2 (f}) is the local minimizer of S ', it is easily to verify that d* satisfies (|3.2p and 
thus is the critical point of E. As a consequence, Corollary 13.11 holds for d* with constants 9, f3 
determined by d*. By ()4.9p . 9' can be determined by 9 and 7. In addition, we further choose 9' 
smaller if necessary such that (recall that 7 > 1) 

9' < X— L (5.28) 
27 

Let us define 

w = mm < — , — , 



2' 2' IOC4 
and set 

t = sup{t G [0,to], ||d(t) -d*|| H i < w, V s G [0,t)} 

If we assume 

cj 2 < -tu, (5.29) 

then by the continuity of d(t) in H 1 ($7), we have to > 0. Next, we shall prove that to > to by 
contradiction. We introduce the auxiliary functional 

+00 



/+00 
r(T)dT, 



and the function 

d(t) = d(t)-d E + d* E . 

It easily follows that 



> E(d(t)) - E(d*) = E(d(t)) - g(d(t)) + g (d(t)) - E(d*) 

= [ F(d(t))-F(d{t))dx + g{d(t))-E(d*). (5.30) 
Jn 
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By definition, d(t) £ J\f. Moreover, on [0,to], 

\\d(t) - d*|| H i < \\d(t) - d*|| H i + \\d E - d* E \\ u i 

5 r +OQ 

< w + c\\h(t)-h 00 \\ H i {r) <- + cJ ||h t (r)|| H 



+00 

1 ._. dr 



3(T) 



5 

< 2 +cM i- 



Taking 



M x < mini 1,^1, (5.31) 

then we have ||d(i) — d*|| H i < 5. Since d* is a local minimizer of S 1 , we see that 

<?(d(t)) - E(d*) = <T(d(*)) - <?(d*) > 0, t e [0, t ]. (5.32) 

On the other hand, since |d(£)| K 3 < 1 and |d(i)| K 3 < 3 (this is due to the maximum principle), 
we infer from the standard elliptic estimate and (H5') that 

r r+00 

/ F{d(t)) - F{d{t))dx < CfeN " d E + d* E \\ < C 5 c / ||h t (r)|| L2(r) dr 

< C b cM^~ l {l + t)~\ (5.33) 

Let us introduce now two further functions 

z(t) = (C 5 c + 1)M 3 7 _1 (1 + t)"\ ¥(t) = *!(*) + 

We deduce from (l53DD - (l533l that 

*(*) > M37 _1 (l + t)- 7 > 0, t€ [0,t ], 

and by the basic energy inequality (|2.16p 

|*(i) = j t £(t)-2r(t)-(C 5 c+l)M 3 (l + t)-^ 

< -lmm{v, l}V 2 (t) - {C 5 c+ 1)M 3 (1 + i)" 1 " 7 



1 



1+- 



< -C 6 [V(t) + Mi(l + t) 2 



2 



where 2?(i) = ||Vv(i)|| + ||Ad(i) - f(d(t))||. Arguing as to get KT2\i . using Remark S3] and 
assumptions (HI'), (H4'), we deduce 

^{t) 1 - 9 ' < C 7 (V(t) + (Mi + M 2 )(l + t)-( 1 - e ')( 1 +T) + Mg 1 "^! + i)-^'^ . 

Assuming 

Ah<^Ml, M 2 <^Ml, M 3 <1, (5.34) 



we can see that 



^(t) 1 " 9 ' < C 7 (v{t) + 2M^{l + t)- {l ~ e '^ 
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As a result, we find 

2 



— 1+T 

C 6 (P(t) + M 3 2 (l + t)- — 



C 7 X>(t) + M 3 2 (1 + i)-(i-^')7 



1+7 



> C 8 \V{t) + M 3 2 (1 + t)~ — j , (5.35) 

where we have used the fact that < (1 - 0')7 (cf. (15381) ). It follows from (ITOD . (I5T23D . 
(Oil) . ([Q5]) . assumptions (HI'), (H4'), (H5') and the definition of ^ that 

l ° \\d t (t)\\dt<C 9 ^(0) e ' 

+00 

r(t)dt + z(0) 







< C10 (j|v || 2 + ||d - d*|| H i + ||d B0 - d||| H i + J 



< C u ^||v || 2 + ||d -d*|| Hl +M 3 2 J . (5.36) 

By (|5.23p . (|5.36p and an interpolation inequality, we get 

||d(t ) - d*|| H i < ||d(t ) - do||Hi + ||d - d*|| H i 
< C 12 (||d(t )|| H2 + ||d || H2 )5||d(t ) - doll* + ||d - d*|| H i 



0' 6' 



Taking now 



< Ci 3 ||v |r +||do-d*|| 2 +M3 4 +||d -d*|| H i. (5.37) 



. . { JT^ ( w \$\ . J 1 f tJ7 

01 < mm < 1, - L — , — — — > , (72 < mm < 1, -ro, 1 



M3 < min < 1, 



6c 13 y j ' [ '4 ' V6C 13 

4 ^ 



we infer from (|5.37p that 

3 

||d(t ) — d*|| H i < -w < w. 
This leads to a contradiction with the definition of to- As a result, we have to > to, and 

||do-d(i)|| H i < ||d -d*|| H i + ||d* -d(t)|| H i <c7 2 + w< -w< Vt G [0,t |. (5.38) 

Thus, (|5.24p holds on [0,io]> which implies 

rto 



(H|Vv(t)r + ||Ad(t)-f(d(t))||^<eo. 







As in Proposition 15.11 there exists t* G [y,io] such that 

^||Vv(^)|| 2 + ||Ad(t,) - f(d(i,))|| 2 < 2e to X . 
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and again we have Ap{t*) < max } K. Taking t* as the initial time for the Cauchy 

problem (|5.13p . we can extend the (unique) strong solution to [0, |io] and its V x H 2 -norm is 
uniformly bounded by the same constant C3 as on [0, to]- Repeating the above argument in 
[0, |to]i we can verify that (|5.24p still holds. By iteration we can show that (|5.24p holds for all 
t > 0. Hence, our conclusion follows from Proposition 15.21 □ 

Finally, we can conclude with the following local stability result: 

Theorem 5.4. Let the assumptions of Theorem 1 5. 3\ hold. Then any global strong solution given 
by Theorem \5.3\ converges in V x H 2 (S7) to a single equilibrium (0,doo) with doo G N <~) H 2 (0) 
such that (^(doo) = <f(d*). In addition, convergence rate estimates similar to the 2D case hold 
provided that g and h fulfill the corresponding hypotheses (i.e., assumptions (HI), (H4), (H5) 
are replaced by (HI'), (H4'), (H5'), respectively). Actually, d* is (locally) Lyapunov stable, and 
in particular, if d* is an isolated local minimizer of £ ' , then it is (locally) asymptotically stable. 

Proof. Arguing as in Section H] we find lim (||v(i)||v + ||d(i) — dooll^) = for some doo G 

t— >t+oo 

M n H 2 (il). Then the estimate on the convergence rates can be obtained following the proof 
of Theorem 14.21 Recalling the proof of Theorem 15. 3( we actually showed that, for all t > 0, 
||d(i) — d*|| H i < w < tj/3, which implies that (let t be large) 

II doo - d*|| H i < ||d(t) - doollni + ||d(t) - d*|| H i < 2zu. 

Thus, taking d = d^ and if; = d* in Corollary 13. 11 we see that \S'(d 00 ) — £{d*)\ l ~ e = ^(doo) — 
E{d*)\ l ~ e < || - Ad* +f(d*)|| = 0. Therefore, doo is also a minimizer of S because from 
the definition of w, ||doo — d*|| H i < 5. Actually, the proof of Theorem 15.31 implies that d* is 
(locally) Lyapunov stable. Moreover, it is easy to see that if d* is an isolated local minimizer 
then doo = d*. □ 



6 Appendix 

We report some properties of the lifting functions d^ and dp (cf. (j2.7j) and (|2.2ip ) that have 
been used in the previous sections. Below we denote by c a generic positive constant which 
depends on n and f2 at most. 

Lemma 6.1. For any t > 0, and k = 0, 1, 2, j = 0, 1, we have 

(^ \\d{d E m^ <4diHt)\\ nk -^ ry 

(ii) ||dp(t) — d*|| H fc < c][h(t) — h*|| fe _i , where d* is the unique solution to 

f-Ad, = 0, xeQ, 
[d* = h*, x G T. 

Proof. The conclusion follows from the classical elliptic regularity theory (cf., e.g., |20l l28j ) . □ 

Lemma 6.2. Let do G H 2 (fi) with |do|]R« < 1- Suppose that h satisfy ([2.4|) - (|2.5|) and \\t G 

I?, L([0, +00); H2 (r)). Then, for any t > 0, the following estimates hold 



\d P (t)-d E (t)\\h < ce~* [ e T \\h t (r)\\ 2 , dr, (6.2) 

Jo H 



32 



\\d t d P (t)f + \\d P (t) - d E (t)f u2 < c f \\h t (r)\\l h dr, (6.3) 

fwVAdpfdr < c f \\h t (r)\\ 2 h dr. (6.4) 
Jn Jo w 2 ^] 



In addition, we have 

(i) ifh t G L 2 (0,+oo;H^(r)) then 



lim ||d t d P (i)||=0, (6.5) 

t— >+oo 



(ii) ifh t satisfies (H6) then, for all t > 0, 

\\d t dp(t)\\ 2 + \\d P (t)-d E (t)\\h < c {l + t)- 2 - 2 \ (6.6) 

/ \\VAd P (r)\\ 2 dr < c(l + t) 
"'I 

Proof. It follows from flZ2D and (|2~2T|) that 
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(6.7) 



-A(dp - dp) = -dtdp, in x K+, 
d P -dp = 0, onrxl+ 



(6.8) 



and 

'dt{&P - dp) - A(dp - dp) = -d t d E , in n x R+, 
d P -dp = 0, onrxE+, (6.9) 
^dp - dp| t=0 = 0, in U. 

Multiplying the first equation in (|6.9p by (dp — dp) — A (dp — dp), integrating by parts and 
using the Poincare inequality, we obtain 

\j t (\\d P - dp|| 2 + ||V(d P - dp)|| 2 ) + ||V(dp - dp)|| 2 + ||A(d P - dp)|| 2 

< (||dp-dp|| + ||A(dp-dp)||)||a t dp|| 

< (Cp||v(dp - dp)|| + ||A(dp - dp)||)||a t dp|| 

< i||V(d P - dp)|| 2 + i||A(dp - dp)|| 2 + Qc 2 + ||^dp|| 2 , (6.10) 
which, together with Lemma 16. 11 implies 

\\d P (t) -d B (t)\\h < «^ Cli t e c n\dtd E (T)\\ 2 dr < ce- Cl * f e^\\h t (r)\\ 2 , dr, (6.11) 

jo jo H 2 U) 

that is, 

Applying now the Laplacian to the first equation in (|6.9p . we get 



d t A(d P - dp) - A 2 (dp - dp) = 0, in x K+, 

A(dp-dp) = h t , onrxl+, (6.12) 
A(d P -dp)|t =0 = 0, in a 

Multiplying the first equation of (|6.12|) by A (dp — dp) and integrating by parts, we get 
^||A(dp-dp)|| 2 + ||VA(dp-dp)|| 2 
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< ||a n A(dp-d E )|| H _, (r) ||h t || H x (r) < c ||A(dp-d B )|| Hl ||h t || H i (r) 

< i(||VA(d P - d E )|| 2 + ||A(d P - dp)|| 2 ) + c\M 2 Hh(r y (6-13) 
Hence, from (|6.10)) and (|6.13j) we infer 

j t \\d P (t) - d E (t)\\h + c 2 (\\d P (t) - d E (t)\\lv + l|VA(dp - d E )\\ 2 ) < cllhtl^ (6.14) 



which entails (16.41) and 



\\d P (t)-d E (t)\\ 2 H2 <c [ \\h t (r)\t , dr. 

JO ^2(1) 



(6.15) 



Thus (Q follows from ffTTBT) and the fact ||$dp(*)|| = ||Ad P 
Now if hi G L 2 (0, +oo; Ha(r)), we infer from (lOOl) that 

/ ||A(d P (t) - d E (t))\\ 2 dt < c / \\d t d E (t)\\ 2 dt < c / ||h t (t)|| 2 i dt < +oo. (6.16) 
Jo Jo Jo H ^(r) 

Then it follows from (IBTTSD . (I6J61) and Lemma O that lim^ +00 ||A(d P (t) - d E (t))\\ 2 = 0, 
which implies (|6.5[) . 

Furthermore, if (H6) holds, then (|6.14p implies that (cf., e.g., 



||d P (t) - d E (t)\\ 2 H2 < c (i + ty 2 ' 2 \ v t > o. 

Using (|6.14p once more, we have 

T \\VAd P (T)\\ 2 dr = I ||VA(d P - d E )( T )\\ 2 dT 

•* 2 ^ 2 



< c 



n -2-2 7 c , ^-l-2 7 



< c^i + -j + i+^i+ 5 J ^ ( i +r ^, vt> , 

and this gives ()6.7|) . The proof is complete. □ 
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